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Summary 

It  is  shown- that,  by  a  reparametrization,  the  problem  of  estimating  a  linear 
combination  of  variance  components  can  be  reduced  to  that  of  estimating  a  single 
variance  component.  Such  a  reduction  is  used  to  obtain  some  characterizations 
of  nonnegatively  estimable  linear  combinations  of  variance  components.  Charact¬ 
erization  of  nonnegative  estimability  using  MINQUE  is  also  discussed. 
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Key  words  and  phrases.  Nonnegative  estimability.  MINQUE  (given  Z) ,  MINQUE  (E,NND) 

quadratic  subspace,  subspaces  that  preserve  nonnegative 
definiteness . 


1.  Introduction.  Recently,  much  attention  has  been  given  to  the  problem  of 
nonnegative  estimation  of  variance  components.  Several  nonnegative  estimators 
have  been  proposed  sacrificing  unbiasedness;  see  e.g.  ,  P.S.R.S.  Rao  and  Chaubey 
(1978),  Hartung  (1981)  and  Chaubey  (1983).  In  the  works  of  LaMotte  (1973), 
Pukelsheim  (1981a, b),  Mathew  (1984)  and  Baksalary  and  Molinska  (1984),  the  main 
concern  is  the  existence  of  nonnegative  definite  (nnd)  quadratic  estimators  that 
are  also  unbiased.  The  present  work  is  concerned  with  the  existence  and  character¬ 
ization  of  nnd  quadratic  unbiased  estimators  (QUE's)  of  a  linear  combination  of 
variance  covariance  components . 

The  problem  of  nonnegative  estimation  of  variance  components  is  not  fully 
resolved  by  characterizing  linear  combinations  of  variance  components  that  admit 
nnd  QUE's.  It  should  also  be  possible  to  obtain  an  nnd  QUE  having  some  optimal 
properties  (for  e.g.  properties  similar  to  those  of  C.R.  Rao’s  MINQUE) .  The  work 
of  Pukelsheim  (1981a)  is  a  significant  achievement  in  this  regard.  Under  a 
quadratic  subspace  condition,  Pukelsheim  showed  that  in  order  to  verify  non¬ 
negative  estimability ,  it  is  enough  to  check  the  nonnegativity  of  the  MINQUE  (given 
I) .  This  actually  solves  the  problem  of  nonnegative  estimability  of  variance  com¬ 
ponents  from  balanced  data,  since,  as  observed  by  Anderson,  et.  al.  (1984  p.  170), 
the  quadratic  subspace  condition  is  always  satisfied  in  this  case.  Pukelsheim's 
result  has  been  extended  by  Mathew  (1984).  However,  the  procedure  outlined  in 
Pukelsheim  (1981a)  and  Mathew  (1984)  will  not  always  work  for  unbalanced  data. 

In  the  next  section  we  show  that,  by  a  suitable  reparametrization,  the 
problem  of  estimating  a  linear  combination  of  variance  components  can  be  reduced 
to  that  of  estimating  a  single  variance  component.  This  reduction  has  enabled 
us  to  obtain  some  characterizations  of  nonnegatively  estimable  linear  combinations 
of  variance  components  and  also  to  obtain  nnd  QUE's,  A  complete  solution  is 
given  to  the  nonnegative  estimation  problem  in  the  case  of  a  model  with  two 


variance  components.  These  are  discussed  in  section  3.  In  section  4  we  consider 
the  problem  of  characterizing  nonnegative  estimability  using  MINQUE. 


2.  Notations  and  Preliminary  Results. 

Let  Y  be  a  random  Rn-vector  with  E(Y)=  X8  and  D(Y)=  V0 


i=l 


eiV 


Here  X 


is  a  known  n*m  (m<n)  matrix,  8  is  a  vector  of  unknown  parameters  varying  over 
Rm,  V^(i=  1,2,. . . ,k)  are  known  real  symmetric  matrices  and  9=  (9^ ,0^ , . • • ,9^) ' 
is  a  vector  of  unknown  parameters  varying  over  the  set  0,  a  subset  of  R  . 

The  following  assumptions  are  made  regarding  the  matrix  Vfl  and  the  set  0 


(i) 

for 

each  0  6  0,  V0 

is  nnd 

(1) 

(ii) 

the 

elements  of  0 

_k 

span  R 

(2) 

(iii)  there  exists  an  nnd  matrix  V  €  SP{V  :  0 e  0}  such  that 

o  u 

for  i-  1,2, . . . ,k  R(V±) c  R(Vq) 

Here  R(*)  denotes  range.  We  denote  the  above  model  as 


(3) 


Y~  (X8,  V0),  06  0  (I) 

The  unknown  0^'s  could  be  components  of  variance  or  covariance.  We  are  interested 
in  estimating  a  linear  combination  q’9  =  2+* '  ’+<^k®k'  tlle  estimators  under 

consideration  heing  quadratic  forms  in  Y.  We  assume  without  loss  of  generality 
that  q’q  =  1  and  q^O.  For  an  n>n  positive  definite  (p.d.)  matrix  Z,  let  M^.  = 

I-  X(X'Z  ^X)  X’Z  ^  (A  denotes  a  generalized  inverse  of  the  matrix  A)  and  let 
Sj.  be  the  matrix  whose  ( i j ) element  is  tr  Z  (i,j*  l,2,...,k). 

Then  it  can  be  shown  that  R(Sj.)  does  not  depend  on  Z  and  q‘9  has  an  invariant 
quadratic  unibased  estimator  (IQUV)  iff  qe  R(S^) .  When  qe  R(S^.),  the  MINQUE 
(given  Z)  of  q'9  is  given  by  Y'(  J  A  Z’m^V  M^Z"  )Y,  where  A-  (A^,...^)'  is 


any  solution  to  S^.X=  q.  The  MINQUE  (given  £)  of  q*0  is  the  unique  estimator  ob¬ 
tained  by  minimizing  tr  AZAT,  where  A  is  any  symmetric  matrix  satisfying  the  con¬ 
ditions  AX=  0  and  tr  AV^=  q^  (i=  1,2, . . . ,k) .  For  the  details,  we  refer  to 
Rao  (1970,  1971,  1972,  1973),  Kleffe  (1977)  or  Rao  and  Kleffe  (1980).  For  our 
discussion,  it  is  conveneint  to  have  the  following  definitions.  Definition  1 
is  due  to  Pukelsheim  (1981a). 

Definition  1  We  say  that  q'9  is  nonnegatively  estimable  in  model  (I)  if  it  has 
a  nonnegative  definite  quadratic  unbiased  estimator. 

Definition  2  We  say  that  MINQUE  (given  £)  characterizes  nonnegative  estimability 
in  model  (I)  if,  for  every  vector  qe  R  ,  the  nonnegative  estimability  of  q'0 
implies  the  nonnegativity  of  its  MINQUE  (given  £) . 

For  a  nonnegatively  estimable  q'0  and  for  a  p.d.  E,  let  minimize  tr  BZBZ, 
where  B  is  any  nnd  matrix  such  that  Y'BY  is  an  unbiased  estimator  of  q'0.  We 
shall  refer  to  Y'B*Y  as  the  MINQUE  (I,  NND)  of  q'0.  Let 


Q  = 


'  2 

i-91  -qLq2  -qxq3  . "Vk-1  ^ 

— q2qi  ’VS  . -q2qk-l  q2 


~qk-lql  "qk-lq2  _qk-iq3 .  1  qk-l  qk-l 

"qkql  "qkq2  -qkq3  .  -qkqk-l  qk 


Then  |q|  =  q^,  which  is  nonzero  (by  our  assumptions  q’q=  1  and  q^  0).  For 
n  =  (ni,n2,...,nk),t  if  we  let  0  =  Qn,  then  e1v1+e2V2^' • ,+0kVk=  ni^Vl_qlVq)  + 

vw^-n-A-A-iVnv  Here  vq  ■  J1qivi- 


We  now 


consider  the  model 


I 


4 

k-1 

Y~  (xe,  _Z  ni(V1-qiVq)4fikVq)  (II) 

where  r\e  Q  10  =  { Q~ 1© :  9  e  0}. 

Lenana  1  (i)  Every  QUE  (or  IQUE)  of  q'0  in  model  (I)  is  a  QUE  (respectively  IQUE) 

of  nk  in  model  (II)  and  vice  versa. 

(ii)  The  nonnegative  estimability  of  q'0  in  (I)  is  equivalent  to  the 
nonnegative  estimability  of  nk  in  (II). 

(iii)  For  any  p.d.  E,  the  MINQUE  (given  E)  of  q'0  in  (I)  is  same  as  the 
MINQUE  (given  E)  of  nk  in  (II). 

(iv)  If  q'0  is  nonnegatively  estimable  in  (I),  then  for  any  p.d.  E,  the 

MINQUE  (E,  NND)  of  q'0  in  (I)  is  same  as  the  MINQUE  (E,  NND)  of  r)^  in  (II). 

Proof:  Y'AY  is  a  QUE  (or  IQUE)  of  q'0  in  (I)  iff  tr  AW±  =  qt  (i=  1,2 . k) 

2 

and  X'AX  =  0  (respectively  AX  =  0) .  Then  tr  AV  =  Eq.  =  1  (by  assumption). 

1  k-1  q 

Since  V,  -q. V  * - T  q.(V  -q.V  ),  the  condition  tr  AV  *  q  (i=  1,2,. . . ,k) 

k  k  q  qk  l  i  i  q  i  i 

is  seen  to  be  equivalent  to  the  conditions  tr  AV  =1  and  tr  A(V.  -  q.V  )  =  0 

q  l  l  q 

(i =  1,2 , . . . ,k-l) .  This  proves  the  assertion  in  (i) .  (ii)  follows  from  (i). 

Since  the  class  of  IQUE's  of  q'0  in  (I)  and  in  (II)  are  the  same,  the  minimum 
norm  element  in  this  class  gives  the  MINQUE  (given  E)  of  q'0  as  well  as  n^.  This 
proves  part  (iii).  (iv)  follows  similarly.  □ 

3.  Nonnegative  Estimation. 

When  the  matrices  V^  (i=  l,2,...,k)  are  nnd,  conditions  for  the  nonnega¬ 
tive  estimability  of  a  single  variance  component  in  the  model  (I)  has  been  derived 
by  Pukelsheim  (1977,  Theorem  5.1).  Altemativef orms  of  the  same  condition  are 
given  in  Kleffe  (.1977,  Theorem  3)  and  Rao  and  Kleffe  (1980,  Theorem  5.5.1).  We 
now  proceed  to  obtain  similar  results  for  the  nonnegative  estimability  of  q'0. 


The  following  lemma  will  be  used  in  the  sequel.  When  a  symmetric  matrix  A  is 
nnd,  we  denote  A>0. 


Lemma  2  (Pukelsheim,  1981a).  Let  M=  I-XX  .  Then  q'0  is  nonnegatively  estimable 
in  (I)  iff  t^MVjM+. .  .+t^MV^M  0  implies  q^t^+. .  .+q^  >  0 . 

Let  P^  denote  the  orthogonal  projector  onto  the  subspace  R(HV^M-q^MV  M)+... 
+R(MV^  ^M-q^  ^MV  M) ,  i.e.  if  the  columns  of  the  matrix  H  form  a  basis  for  this  sub¬ 
space,  then  P^=H(H'H)  H' .  Further,  let  l/^(q)  denote  the  vector  space  spanned 
by  the  matrices  MCV^q^V  )M  (i=  1,2,.  . .  ,k-l) . 

Theorem  1.  (i)  If  (I-P^)MV^M  (I-P^)  is  non-null,  thenq'0  is  nonnegatively  estimable 

in  (I)  iff  (I-P  )MV  M(I->P  )  is  nnd. 

q  q  q 

(ii)  If  there  exists  an  nnd  matrix  WQ e  i/^q)  satisfying  R(MV^M-q  JfV^M) c  R(WQ) 
for  i=  l,2,.,.,k-l,  then  q'0  is  nonnegatively  estimable  in  (I)  iff  (I-P  )  MV  M(I-P  ) 

q  q  q 

is  nonnull  and  nnd. 

Proof :  Let 


U,  -  v,-q4V  for  i  =  1,2 ... .  ,k-l 
i  i  i  q  ’ 

=  V  for  i  =  k. 

q 

From  Lemma  2  and  Lemma  1  (ii),  it  follows  that  q'9  is  nonnegatively  estimable  in 

(1)  iff  t ^MU ^M+ .  . .  +t^MU^M  >_  0  =>  t^_>  0,  where  M=  I-XX+.  In  view  of  assumption 

(2)  about  0,  there  exists  a  nonnull  t,  satisfying  t  MU  M+. . ,+t  MU  M  >  0.  This 

k  11  k  k.  — 

implies  t.  (I-P  )  MV  M(I-P  )>  0.  If  (I-P  )MV  M(I-P  )  is  a  nonnull  matrix, 

v  k  q  q  q  —  q  q  q 

since  t.  is  nonnull,  t.  is  positive  iff  (I-P  )MV  M(I-P  )  is  nnd.  This  proves 

iC  K.  q  q  q 

part  (i).  The  "if"  part  in  (ii)  is  clear  from  (i) .  The  "only  if"  part  also 
follows  from  (i;  once  it  is  shown  that  for  q'0  to  be  nonnegatively  estimable, 
(I-P  )MV  M(I-P  )  must  be  nonnull  when  there  exists  a  matrix  W„  as  specified 

q  q  q  v  0 


Suppose  ( 1-P  )MV  M(r-P  )  =  0.  Let  J  t.MU.M  >  0  where 

q  q  q  i^1  i  1 


in  the  theorem. 


6 


t  ^  0.  Since  Pq  MU.M=MU.M  for  i=  l,2,...,k-l,  F  t.MU.M^O  can  equivalently  be 

k-1  1  1  i=l  1  1 

written  as  P  (  F  MU.M)P  +  t.  [ (I-P  )MV  MP  +  P  MV  M(I-P  )  +  P  MV  MP  ]  >0.  This 
wrrccen  as  r  k  L  i  q  kl  q  q  q  q  q  q  q  q  q  - 

1=1 

gives  (I-P  )MV^MP^=0,  which,  along  with  (I-P^MV^MCI-P^)  =  0  yields  (I-P^)MV^M=0 

or  equivalently  R(MV^M)  c  R(Wq)  and  hence  there  exists  t^  <0  such  that 

UL+t,  MV  M  >  0.  This  contradicts  the  nonnegative  estimability  of  q'0  since 
0  k  q  — 

Wq  is  a  linear  combination  of  MU^M  (i=  1,2, . . . ,k-l) .  U 

Remark  1  If  (I-P  ) MV  M(I-P  )  is  nonnull,  then  it  cannot  be  indefinite.  From  the 
-  q  q  q 

proof  of  Theorem  1  it  follows  that  there  exists  t^?4  0  satisfying  t^(I-P  ) MV^M(I-P^) 
>  0.  Hence  (I-P  )  MV  M(I-P  )  is  either  nnd  or  nonpositive  definite.  In  case  it 

-  q  q  q 

is  nonpostive  definite,  q'0  has  a  nonpositive  definite  quadratic  unbiased  esti¬ 
mator. 

In  the  case  of  a  model  with  two  variance  components,  the  following  corollary 

gives  a  complete  characterization  of  nonnegatively  estimable  q'0. 

Corollary  1.  Suppose  k=  2  in  model  (I).  Let  P  and  M  be  as  defined  before. 

-  q 

(i)  If  (I-Pq)  MV^M(I-P^)  is  nonnull,  then  q'0  is  nonnegatively  estimable  iff 
(I-P  )  MV  M( I-P  )  is  nnd. 

q  q  q 

(ii)  If  MV_^M-q^MV^M  is  nnd  or  nonpositive  definite,  then  q'0  is  nonnegatively 
estimable  iff  (I-P  )  MV  M(I-P  )  is  nonnull  and  nnd. 

q  q  q 

(iii)  If  MV  M-q.MV  M  is  indefinite  and  if  (I-P  )  MV  M(I-P  )=0,  then  q'0  is 

1  1  q  q  q  q 

nonnegatively  estimable  iff  (a)  R(MV^M) c  R(MV^M-q^MV^M)  and  (b)  there  exists  a 

real  number  a  such  that  MV  M+a(MV,M-q  MV  M)  is  nnd. 

q  1  1  q 

Proof :  Only  part  (iii)  needs  to  be  proved.  Suppose  q'0  is  nonnegatively  esti¬ 


mable.  If  (I-P  )  MV  M(I-P  ) =  0,  then  proceeding  as  in  the  proof  of  Theorem  1  (ii), 

q  q  q 

we  get  R(MV  M)  c  R(MV1M-q.MV  M) .  Since  there  exists  a  nonnull  real 

q  1  L  q 


number  satisfying  t^ (MV^M-q^MV^M) +  t2MV^M  >_  0 ,  the  nonnegative  estimability  of 
q'9  demands  that  there  exists  a  positive  t£  satisfying  the  same,  which  leads  to 
part  (b) .  To  prove  the  sufficiency  of  the  conditions,  we  observe  that  since 


MV^M-qfMVqM  is  indefinite,  conditions  Cal  and  (b)  guarantee  the  existence  of  an 

nnd  matrix  A  satisfying  tr  A(MV_M-q  MV  Ml=  0  and  tr  A(MV  M+a(MV,M-q.,MV  M))  >  0 

llq  qllq 

or  equivalently  tr  AMV  M  >  0.  Then  - ~^r— rr  Y'MAMY  is  an  nnd  unbiased  estimator 

7  q  tr  AMV  M 

q 

of  q'0.  □ 


In  the  case  of  a  model  with  two  variance  components  with  V^,V2  nnd>  we  shall 
now  obtain  explicit  characterization  of  vectors  q  for  which  q'9  is  nonnegatively 
estimable.  To  this  end,  let  X^>_  '  •  *  jl  denote  the  proper  eigenvalues  of 

MV^M  w.r.t.  MV^M  (see  Rao  and  Mitra  (JL9.71,  Section  6,3),  Here  r  is  the  rank  of 
MV^M.  Some  of  the  X^'s  could  be  zero  and  some  of  them  may  be  repeated. 

Corollary  2  Let  k  =  2  in  model  (I)  and  suppose  and  are  nnd.  Let  X^(i=  1,2,... 
be  as  defined  above. 


(a)  Suppose  R(MV^M) n  RCMV^M)  =  {0 } .  If  MV^M  and  MV^M  are  nonnull,  then  q'9  is 
nonnegatively  estimable  iff  q^>0,  q2>0. 

(b)  Suppose  R(MV1M)  n  R(MV2M)  j*’  (0).  Then 

(i)  if  R(MV^M)  i  RfMV^M) ,  q'0  is  nonnegatively  estimable  iff  q  _>q2Xr_>0. 

(ii)  If  R(MV^M) c  RCMV^M)  and  rank(MV^M)  <rank(MV2M),  then  q'0  is  nonnegatively 
estimable  iff  q2X^^q^>0. 

(iii)  If  R(MV^M) =  R(MV2M),  then  q'0  is  nonnegatively  estimable  iff 

Vl-ql-q2Ar-°- 


Proof:  Since  and  V 2  are  nnd,  any  nonnegatively  estimable  q'0  must  have  q^j>0, 
q2>0  (LaMotte,  1973).  Using  Theorem  6.3.4  in  Rao  and  Mitra  (1971),  we  see  that 
there  exists  a  nonsingular  matrix  P  satisfying  P'MV2MP=  diag  (1^,0)  and  P'MV^MP= 
diag(A1,A2)  where  A^  diag(X  1*X2,  •  •  ■  »^r)  and  A2  =  diagfX^, . . .  ,*n) . 
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(a)  R(MV^M)  n  =  {0}  iff  A^=0.  In  such  a  case,  for  q^>0,  q2_>0, 

ql  +  q2 

rank (\ — J  2 ^  +  ~  ^ ^  diag(I^_  >0)P  'Y  is  an  nnd  unbiased  estimator  of 

q'0. 

(b)  MVqM=  P'  1diag(X1q1+q2>.  ..  ,Arq1+q2,q1Ar+1...  .,q1An)  and  MV^-q^MV^  q2P '  1 

2  2 

diag(q2A1~q1, . .. ,q2Ar-q1>q2Ar+1, . . . ,q2A  )  (using  the  assumption  q. +q2 =  1) .  If 
or  q2  is  zero,  then  it  can  be  verified  that  the  necessary  and  sufficient  con¬ 
dition  for  the  nonnegative  estimability  of  a  single  variance  component  as  given  in 
Theorem  5.1  of  Pukelsheim  (1977)  or  Theorem  3  in  Kleffe  (1977)  or  Theorem  5.5.1 
in  Rao  and  Kleffe  (1980)  reduce  to  those  given  in  the  theorem.  We  now  consider 
the  case  f  0,  q2  ^  0. 

(i)  When  R(MV.M)  4-  R(MV  M) ,  A.  4  0.  Let  us  assume  A  >  0,  If  q  A  -  q  >  0,  then 
1  22  r+1  2  r  1 

MVlM”qlMVqM  nnC*  and  ^as  t^le  same  range  as  MV^M.  Hence  from  Corollary  l(ii) 

it  follows  that  q'0  is  not  nonnegatively  estimable.  Conversely,  if  q^^-q^  — 

since  A  , .  >  0  and  MV  M>  0,  we  can  find  an  nnd  A  satisfyingtr  AMV  M>  0  and 
r+1  q  —  9 

tr  A(MV  M-q  MV  M)  =  0.  Then  - „  Y'MAMY  is  an  nnd  unbiased  estimator  of  q'8. 

1  1  q  tr  AMV  M 

q 

We  observe  that  R(MV^M)  c  R(MV2M)  iff  A2=0  and  then  rank  (MV-^M)  <rank(MV2M)  iff 
A^  =  0.  Using  these  observations,  the  rest  of  the  corollary  can  be  established 
using  Corollary  1.  □ 

n_.  Corollary  2  reduces  to  Theorem  1  in  Baksalary  and  Molinska  (1984)  when 


V 


2 


I. 

Let  Vq  be  the  nnd  matrix  as  defined  in  assumption  3.  Let  Uq  be  an  nnd  matrix 


satisfying  MVqMUq  =  0  and  MVqM  +  Uq  is  p.d. 

Theorem  2.  Suppose  (I-P  )MV  M(I-P  )  is  a  nonnull  nnd  matrix  of  rank  rn .  Then 

-  q  q  q  0 

(i)  the  estimator  —  Y'[(I-P  )MV  M(I-P  ))  Y  Is  an  nnd  unbiased  estimator  of  q'0. 

q  q  q 

(ii)  if  there  exists  an  nnd  matrix  Wq  as  specified  in  Theorem  1  (ii)  ,  the  above  es¬ 
timator  is  the  MINQUE(MVqM+Uq ,  NND)  of  q'0. 


Proof:  (i)  We  observe  that  tr[((I-P_)MV  M(I-P_)]  (V,-q,V  )*0  for  i*  1,2,..., k-1 

-  q  q  q  i  i  q 

and  tr[ (I-P  )MV  M(I-P  ) ]+V  -  tr[(I-P  )MV  M(I-P  )]+[(I-P  )MV  M(I-P  )]  - 

qqqq  qqq  qqq 

rank[(I-P  )MV  M(I-P  )].  Thus,  the  estimator  given  in  the  theorem  is  an  nnd  unbiased 

qqq 

estimator  of  in  (II)  and  hence  that  of  q’0  in  (I). 

(ii)  Let  Y'B^Y  denote  the  nnd  estimator  of  q'0  given  in  part  (i)  of  the  theorem  and 

let  Y'BY  be  any  other  nnd  unbiased  estimator  of  q'0.  We  first  show  that 

B  *  (I-P  )MBM(I-P  )  when  there  exists  a  matrix  Wrt  as  specified  in  Theorem  1  (ii). 
q  q  U 

The  matrix  B  satisfies  BX  *  0,  tr  BV  =1  and  tr  B(V  -q.V  )  *  0  (i  ■  1,2 . k-1) . 

q  i  i  q 

Since  BX=  0  iff  B  =  MBM,  the  last  set  of  conditions  give  tr  B(MViM-q^MV^M)  *  0 

(i  *  1,2, . . .  ,k-l) ,  which  implies  tr  BW^-O  or  equivalently  BW^-  0.  Thus  B 

satisfies  BX-  0,  B(MVM-qJMV  M)  »  0  and  hence  B  -  (I-P  )MBM(I-P  ).  Now  let 

||  B || g  =  tr  B(MV0M-«J0)B(MV0M-KJ0)  -  ||  hJ2Q  +  ||  B-B  J  J  +2  tr  B*(MVQM-HJ0)  (B-B*) 

(MVqM+U0)  -  ||  Bjg+  ||  B-Bj|g  +  2  tr  B*MVQM(B-B*)MV0M.  The  proof  is  complete  if 

we  show  that  the  last  term  is  zero.  Since  V*e  sp{V  ,V  , . . . ,V  }*  sp{V  ,V  -q  V . 

2.  ^  K  qxiq 

V,  -q.  V  },  V.-a.V  +  7  a.(V  -q.V  ),  where  a.  0  (since  (I-P  )MV  M(I-P  )J‘0). 

k— 1  k-1  q  0  1  q  ii  iq  1  qq  q 

Then  tr  B.MV  M(B-B.)MV  M  -  a?(trB.MV  M(B-B.)MV  M  and  trB.MV  MBMV  M  -  trB. (I-P  )MV  M 
*0  *  0  1  *  q  *  q  *  q  q  *  q  q 

(I-P  )MBM(  I-P  )MV  M(  I-P  )=—  tr  (I-P  )MV  M(I-P  )MBM  -  —  tr  BV  -  — -  tr  B.MV  MB  MV  M.  □ 

q  qqqrQ  qqq  rQ  q  r0  *  q  *  q 

Remark  3 .  If  we  are  interested  in  the  nonnegative  estimability  of  a  single  var¬ 
iance  component,  Theorem  2  (ii)  reduces  to  Theorem  5.5.2  in  Rao  and  Kleffe  (1980) 
when  the  matrices  V±  (i=  l,2,...,k)  in  (I)  are  nnd. 

2 

As  an  example,  we  consider  a  model  with  two  variance  components  Y  ~  (XB,  0^ 

2  2  2 
Ifa  ®  ^k^k  +  °2  *b  8  ^k^  where  ^  is  a  k-component  vector  of  ones,  o^0,  c^^O.- 

Such  a  model  arises  in  the  interblock  analysis  of  a  block  design  with  b  blocks, 

each  block  having  k  plots.  A  model  from  genetics  discussed  by  Gnot  and  Kleffe 

(1983,  p.  275)  is  also  of  the  above  form.  An  explicit  characterization  of  q^ 


I 


2  2 

and  such  that  q^O^  +  ^2^2  *"s  nonne8at^ve^y  estimable  can  be  obtained  using 

Corollary  2  or  Theorem  1  in  Baksalary  and  Molinska  (1984) .  We  shall  consider  the 

2  2  2 

nonnegative  estimation  of  0^  and  (ka^’^).  These  are  the  parametric  functions 

of  interest  when  we  want  to  combine  the  inter  and  intra  block  estimates  of  the 

2 

treatment  contrasts  in  a  block  design.  For  the  estimation  of  0^,  Theorem  1  (ii) 

and  Theorem  2  (ii)  clearly  applies.  Let  M  =  I  -  XX*,  V,  ■  L  fl  1,  1’  V  =  I,  SI., 

1  b  k  k  l  b  k 


-  and  q„  =  -  .  Then  V  =  - (kV,+V.)  and  V  -q,V 

' - 2  2  / - 2  q  / - 2  12  1  1  9 

1+k  1+k  1+k 


■y(V1-kV2)  . 


Since  V^-q^V^  is  nonpositive  definite,  so  is  MV^M-q^MV  M.  Thus  the  nonnegative  es- 
2  2 

timability  of  ka^-K^  can  be  verified  using  Theorem  1  (ii) .  When  it  is  nonnega- 
tively  estimable  (which  is  the  case  for  a  block  design)  its  MINQUE  (I,  NND)  can 
be  computed  using  Theorem  2  (ii).  Explicit  characterization  of  nonnegatively  es¬ 
timable  q ' 0  for  some  ANOVA  models  is  given  in  Pukelsheim  (1979,  1981a). 


4.  Characterization  of  Nonnegative  Estimabilitv  using  MINQUE. 


We  now  proceed  to  obtain  conditions  under  which  MINQUE  (given  E)  characterizes 
nonnegative  estimability  in  model  (I). 

Let  8  be  a  subspace  of  real  symmetric  matrices  of  order  n.  For  an  n><n  posi¬ 
tive  definite  matrix  N,  let  PN  denote  the  orthogonal  projector  onto  8,  where 
orthogonality  is  w.r.t.  the  inner  product  <&,B>  =  tr  ANBN;  A,B  symmetric. 

Definition  3.  We  say  that  8  is  an  N-quadratic  subspace  if  BNBe  8  whenever  Be  8. 
Definition  4.  We  say  that  8  preserves  nonnegative  definiteness  wr.r.t.  N  if 
P,.(B)  >  0  whenever  B>  0. 

Definition  3  is  given  in  Musiela  and  Zmyslony  (1978,  Appendix).  Definition  4 
is  a  generalization  of  a  definition  given  in  Mathew  (1984)  and  is  a  special  case 
of  the  definition  of  a  nonnegativity  preserving  linear  transformation  given  in 
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de  Pillis  (1967).  It  is  clear  that  8  preserves  nonnegative  definiteness  w.r.t  N 
iff  N^2BN^2  =  {N^2BN^2:B e  g}  preserves  nonnegative  definiteness  w.r.t.  I.  Us¬ 
ing  this  observation  characterization  of  subspaces  that  preserve  nonnegative  defi¬ 
niteness  can  be  obtained  similar  to  Lemma  1  and  Lemma  2  in  Mathew  (1984).  It 
can  also  be  shown  that  an  N-quadratic  subspace  preserves  nonnegative  definiteness 
w.r.t.  N,  (cf.  Pukelsheim  1981a,  Lemma  2). 

Let  M^.  be  as  defined  in  section  2  and  let  8^,  denote  the  subspace  spanned  by 
the  matrices  (i  =  l,2,...,k).  The  following  result  is  a  generalization  of 

Theorem  2  in  Pukelsheim  (1981a)  and  the  theorem  in  Mathew  (1984). 

Theorem  3.  (i)  MINQUE  (given  Z)  characterizes  nonnegative  estimability  in  (I)  iff 
8^.  preserves  nonnegative  definiteness  w.r.t.  Z  \ 

(ii)  Let  8^  preserve  nonnegative  definiteness  w.r.t . jT^and  suppose  8^.  is  k-dimension- 


The  theorem  follows  from  the  corresponding  results  for  the  case  1=  I  once 
it  is  observed  that  MINQUE  (given  Z)  characterizes  nonnegative  estimability  in  (I) 
iff  MINQUE  (given  I)  characterizes  nonnegative  estimability  in  the  model 

z"1/2y  ~  (i_1/2xe,z“1/2vez"1/2). 

In  the  introduction,  it  has  been  pointed  out  that  MINQUE  (given  I)  always 
characterizes  nonnegative  estimability  in  a  general  m-way  classification  model 
with  balanced  data.  We  shall  now  apply  the  reuslts  in  this  section  to  the  multi¬ 
variate  linear  model  Y  ~  ((Ip  ®  X)0,  }  8  V),  where  $  of  order  p  is  unknown  and  V  is 
a  known  nnd  matrix.  Let  U  be  an  nnd  matrix  satisfying  (V+XX')U  =  0  and  G=  V+XX'+U 
is  p.d.  Let  M_  *  I-X(X'G  ^X)  X’G  Then  subspace  8„„  consisting  of  matrices 
S8  M_VM'  is  an  I  8  G  1  -quadratic  subspace  of  dimension  P(P~I~^ ,  where  S  is  any 
symmetric  matrix  of  order  p.  Hence  for  checking  the  nonnegative  estimability  of 
a  linear  combination  of  the  components  of  it  is  enough  to  check  the  nonnega¬ 
tivity  of  its  MINQUE  (given  I8G-'*").  Furthermore,  the  estimate  of  $  obtained  from  the 


MINQUE  (given  I  B  G  )  of  its  components  is  nnd.  For  the  case  V  =  I,  these 
observations  are  given  in  Pukelsheim  (1981a,  p.  295). 
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